Schwarzschild radius from Monte Carlo calculation of the Wilson loop 
in supersymmetric matrix quantum mechanics 
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In the string/gauge duality it is important to understand how the space-time geometry is encoded 
in gauge theory observables. We address this issue in the case of the DO-brane system at finite 
temperature T. Based on the duahty, the temporal Wilson loop operator W in gauge theory is 
expected to contain the information of the Schwarzschild radius i?sch of the dual black hole geometry 
as log(W) = -Rsch/(27ra'r). This translates to the power-law behavior log{W> = 1.89- {T/ X^^^)~^/^ , 
where A is the 't Hooft coupling constant. We calculate the Wilson loop on the gauge theory side 
in the strongly coupled regime by performing Monte Carlo simulation of supersymmetric matrix 
quantum mechanics with 16 supercharges. The results reproduce the expected power-law behavior 
up to a constant shift, which is explainable as a' corrections on the gravity side. 
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Introduction. — String/gauge duality, which origi- 
nated from the AdS/CFT correspondence [ij, has been 
investigated intensively over the past decade. Remark- 
able developments that have been achieved include gen- 
eralization to various cases, confirmation by explicit cal- 
culations, and applications to various branches of physics 
such as hadron physics and condensed matter physics. 

From the viewpoint of string theory, the duality en- 
ables us to study quantum aspects of gravity including 
its non-perturbative effects from the gauge theory side, 
which is more tractable. In this regard it is important 
to understand how gauge theory captures the informa- 
tion of space-time geometry 0. Based on the duahty 
at finite temperature 0], one can show that a temporal 
Wilson loop operator in gauge theory is related directly 
to the Schwarzschild radius which is a fundamental 
quantity that characterizes the dual black hole geometry. 
(See also Refs. Q for related works.) In this Letter we 
confirm this prediction by first-principle calculations on 
the gauge theory side. Note that this is the first confirma- 
tion of the prescription for calculating the Wilson loop 
based on the string/gauge duality in a non-conformal the- 
ory without protection by supersymmetry. As such, we 
consider our results to have impact also in applications 
of the duality to realistic gauge theories. 

The string/gauge duality we study is the one asso- 
ciated with a stack of N DO-branes in type IIA super- 
string theory at finite temperature T. The worldvolume 
theory of the DO-branes is given by Id U(-/V) gauge the- 
ory or matrix quantum mechanics (MQM) with 16 su- 
percharges, and the dual geometry is given by the near- 
extremal black 0-brane solution in type IIA supergravity. 
The supersymmetric MQM has been studied by Monte 



Carlo simulation 8] in the Fourier space [9|. In particu- 
lar, the results for the internal energy at various (effec- 
tive) 't Hooft coupling constant interpolated nicely the 
weak coupling behavior obtained by the high tempera- 
ture expansion 10], and the strong coupling behavior 
predicted by the black hole thermodynamics of the dual 
geometry. (Consistent results were obtained also by us- 
ing a lattice approach [lll|.) 

Here we apply this method to the calculation of the 
temporal Wilson loop operator, and demonstrate that 
one can extract the Schwarzschild radius of the dual black 
hole geometry from it. See Ref. ^13] for earlier discussions 
on a similar issue in the same model using other observ- 
ables and other calculation techniques. 

Wilson loop in the dual string theory. — Let us re- 
view the calculation of the Wilson loop based on the 
string/gauge duality 6] for general D-branes. In addi- 
tion to a stack of N D-branes, which are placed on top 
of each other creating a curved background geometry, we 
consider a single probe D-brane, which is placed far away 
from them in parallel. Since D-branes are objects which 
a fundamental string can end on, we may consider such 
a string stretched between the probe D-brane and one 
of the N D-branes. The amplitude for the string prop- 
agating along a certain loop C on the D-brane can be 
calculated in two different ways. 

First in the worldvolume theory of the N D-branes, 
the process is viewed as a heavy test particle in the fun- 
damental representation of the U(A^) group propagating 
along the loop. The amplitude is therefore given by 



A = {W{C)) e 
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where W(C) represents the Wilson loop associated with 
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the loop C, whose perimeter has the length The mass 
M of the test particle, which appears in H]), is given by 
the distance of the N D-branes and the probe D-brane. 

Next we view the same process on the gravity side as a 
string propagating in the curved space-time background, 
which is created by the N D-branes. The amplitude is 
calculated by the path-integral over the worldsheet at- 
tached to the loop C on the probe D-brane as 



(2) 



where 5'string represents the worldsheet action, whose 
bosonic part is given by the Polyakov action 
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(3) 

Here x*^ represents the embedding of the worldsheet into 
the target space with the metric gMN {M, N = 1, . . . 10), 
while i?(2) represents the two-dimensional scalar curva- 
ture defined for the worldsheet metric hab {a,,b = 1,2) . 
The effective string coupling is given as e'^ in terms of 
the dilaton field (j). We have omitted a term in ([3]) de- 
pending on the NS-NS B-field since the background we 
are considering does not have non-zero B field. 

In what follows, we will be mostly interested in the pa- 
rameter region, in which the string coupling is so small 
that we only have to consider the disk amplitude in ([2]). 
We also restrict ourselves to small a' , which corresponds 
to a large string tension, so that the path integral ^ is 
dominated by the saddle-point configuration. In this pa- 
rameter region, one can use the classical solution to the 
supergravity as the background. According to the dic- 
tionary of the string/gauge duality, the parameter region 
corresponds to taking the planar large- iV limit with large 
't Hooft coupling constant on the gauge theory side. 

Equating ([T|) and ^ , we obtain a formula which relates 
the Wilson loop in the strongly coupled gauge theory to 
the string amplitude on the classical background geome- 
try. The explicit check of this formula has been discussed 
only in highly symmetric cases. In the 4d A/" = 4 super 
Yang-Mills theory (SYM), in particular, it is argued that 
the gauge theory computation of a circular Wilson loop, 
which is half BPS, reduces to a matrix integration with a 
Gaussian weig The obtained result indeed 

agrees with the prediction from the gravity side. The 
agreement can be understood also from the scale invari- 
ance of the worldsheet theory near the D3-branes [lB| . 

The DO-brane case. — From now on, let us restrict 
ourselves to the DO-brane case. The gauge theory side is 
described by the supersymmetric MQM 

M f 1 1 

SsQM = y dttr)^-{DtX,f - -[X,,X,]^ 



where Dt = dt — i [A{t)^ ■ ] represents the covariant 
derivative with the gauge field A{t) being an N x N 
Hermitian matrix. The model can be viewed as a Id 
V{N) gauge theory with adjoint matters Xi{t) {i — 
1, • • • ,9) and 'ipa{t) {a — 1, • • • , 16), which are bosonic 
and fermionic matrices, respectively. The extent of the 
Euclidean time direction (3 corresponds to the inverse 
temperature /3 = T~^, and the fermions ipa obey anti- 
periodic boundary conditions. 

We consider the loop C to be winding once around 
the temporal direction. The corresponding Wilson loop 
operator in the gauge theory is given by p| 



W = — tr Pexp 

N 



dt[iA{t) +n,X,{t) 



(5) 



where n is a unit vector in specifying the direction 
in which the probe DO-brane is separated. Note that 
the adjoint scalar appears in ([5]) unlike the definition of 
the Polyakov line since the end of the string is coupled 
not only to the gauge field but also to the adjoint scalar. 
The overall factor 1 /N is introduced to make the quantity 
finite in the planar large- iV limit. We have used the same 
normalization on the right hand side of 

The gravity dual of the supersymmetric MQM is given 
by the near-horizon geometry of the (Euclidean) near- 
extremal black O-brane solution in tyge IIA supergravity. 
In particular, the metric is given by 




(6) 



where f{U) = 1 - U^/U"^ and do = 2'^n^/^r{7/2). The 
Schwarzschild radius and the inverse Hawking tempera- 
ture are given by 
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Let us evaluate the string disk-amplitude ^ in the 
background geometry ([6]). In the a' — > limit, the second 
term in ([3]) can be omitted, and one can replace the string 
action by the Nambu-Goto action 5ng: which is nothing 
but the area of the string worldsheet times the string 
tension. Following the proposal g!|, we consider a string 
worldsheet localized in the direction. Then, due to 
the form of the metric (O, the Nambu-Goto action for 
the minimal area is given by 5ng — -^PiUoo—Uo), where 
Uoo represents the position of the probe DO-brane. 

Since the perimeter i of the Wilson loop in ([1]) is given 
by £ = /3 in the present set-up, we obtain the identity 



\og{W{C)) - f3M 
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Considering that the mass of the test particle on the 
gauge theory side is given by the position of the probe 
DO-brane, it is natural to identify the second terms on 
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both sides of ([8]). This follows also from the prescription 
proposed in Ref. 17| based on T-duality. Thus we obtain 



log(t^(C)) 



f3Uo 
2tt 
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1.89 



T 



-3/5 



(9) 



where we have used ([7|. 

The range of validity. — Let us recall the range of 
validity for the supergravity description By chang- 
ing the target-space coordinates as U = Uqv?/^ and 
t = ^y/doXUg'^^'^T, the metric ^ and the effective string 
coupling e*^ become 



'4 = idy'icf" 



du' 

W) 



dnt 



,4> = (^(rf-iA^)2i/io 
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(10) 
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where f{u) = u^(l 



-14/5 



) . From (fTO|) . one finds 



that the geometry asymptotes at large u to a geo metry 
which is conformally equivalent to AdS2xS* [18[, and 
that the typical length scale of the geometry is given by 
p = (uT/A1/3)-3/10q,/i/2^ rpj^ig gpg^ig should be much 

larger than the string length a'^^'^ for the a' correc- 
tions to the supergravity action to be negligible. Hence, 
uT /\^/^ ^ 1. In this case, the first term in ([3]), which is 
proportional to p^, becomes large, and the semi-classical 
treatment for the string amplitude ([2]) is also justified. 
Note, however, that we have introduced Uoo- Assuming 
that we only need to require C/qo/C^ to be large (but fi- 
nite), we may assume u to be finite as well. Then we 
obtain the condition T/A^/^ < 1. 

We also need to require the effective string coupling e"^ 
to be small. From (|Tl]), we obtain iV-i^/^i ^ T/A^/^ 
noting that m > 1 in our finite temperature set-up. 

a' corrections. — Let us discuss possible subleading 
terms in ([9]) due to a' corrections on the gravity side. 
There are three effects one should consider: (I) the cou- 
pling with the background field represented by the sec- 
ond term in ([3]), (II) a'-corrections to the background 
fields that appear in the action Q , and (III) the quantum 
fluctuation of the string worldsheet including fermionic 
degrees of freedom in evaluating In order to discuss 
the next-leading order terms, we can treat each of these 
effects separately. 

The effect (I) yields a constant term and a logarith- 
mic term with respect to T/X^/'^ in ([9]) as one can see 
from ([3]) and (fTTjl . The constant term includes logA^, 
but this is canceled by the prefactor \/N \n ^ &s \t 
should. The effect (II) can be neglected at this order 

since a'-corrections to the t ype IIA supergravity action 
,/3 



starts only at the a' order [ig|. The effect (III) yields a 
constant term to ( |9l). This effect is discussed also in the 
case of D3-branes [20|. In fact a logarithmic term can 
appear from it as well due to the insertion of the ghost 
zero mode [14]. 



Monte Carlo simulation. — We perform Monte Carlo 
simulation of the model (|4|) and calculate the temporal 
Wilson loop to check the prediction Q. We use the 
Fourier- mode simulation method in which we take 
the static diagonal gauge A{t) = ■gdiag(ai, • • • , un) with 
— TT < tta < TT, and introduce a cutoff A on the Fourier 
modes as X,[t) = E^=-a ^™e*'^"*, where = ^. Su- 
persymmetry at T = 0, which is broken only due to finite 
A, is shown to recover rapidly as A — > cx) in a simpler 
model [!|. The effective 't Hooft coupling constant is 
given by AeS = X/T^. In actual simulation we set A = 1 
without loss of generality, so that high/low T corresponds 
to weak/strong coupling, respectively. 

Integration over the fermionic matrices yields a Pfaf- 
fian Pf Al , which is complex in general. According to the 
standard reweighting method, one uses |PfA^| to gener- 
ate configurations, and includes the effect of the phase 
when one calculates the expectation values. In fact ViJVl 
is almost real positive at sufficiently high T, but the fluc- 
tuation of the phase becomes larger as T decreases, which 
causes the so-called sign problem. It turned out, however, 
that the results of the reweighting method in the temper- 
ature regime where the sign problem is not so severe are 
actually in good agreement with what we obtain by sim- 
ply neglecting the phase. We interpret this as an effect 
of the large- A^ limit, in which the fluctuations of single 
trace observables vanish. For the same reason, it is ex- 
pected that log(|VF|) agrees with (log |W^|) in the large- A^ 
limit. We therefore calculate the latter in an ensemble 
generated with |PfA^|. Complete justiflcation of these 
simpliflcations is left for future investigations. 

We evaluate ([5]) as a limit W = hm^_»oo W,,, where 



k=G 



M^. = ^tr[] l + ^{iA + n,X,{tk)] 



(12) 



with tk — ^P. The matrices Xi{ti;) are obtained as the 
inverse Fourier transform of the conflgurations generated 
by our simulation. Using the asymptotic behavior W^, ~ 
cogst. large i', we can make a reliable extrapolation 
to J/ = CX3. As the unit vector n, we have used the ones 
in all 9 directions with plus or minus sign in front, and 
averaged over them to increase statistics. 

In Fig. □ we plot (loglM^I) against y-^/s. As T de- 
creases (to the right on the horizontal axis), the data 
show a clear linear growth with a slope consistent with 
the value 1.89 predicted in ([9]). In fact we can fit our data 
to (log \W\) = 1.89T-3/5 _ C, where C = 4.95 for = 4 
and C = 4.58 for A^ = 6. The data points for A^ = 8 are 
very close to those for A^ = 6. Note that the constant 
term and the logarithmic term predicted from the gravity 
side are difficult to distinguish numerically. We therefore 
consider that the value of C extracted above actually 
represents the sum of the two terms at the temperature 
regime investigated. 
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FIG. 1: The plot of (log |W^|) for A = 1 against The 
cutoff A is chosen as follows; A = 12 for TV = 4; A = 0.6/r for 
iV = 6, 8; A = 4 for AT = 14; A = 6 for iV = 17. The dashed 
line represents the results of the high-temperature expansion 
up to the next-leading order with extrapolations to A'^ = oo, 
which are obtained by applying the method in Ref. [l3| ■ The 
solid line and the dotted line represent fits for N — 6 and A'^ — 
4 respectively, to straight lines with the slope 1.89 predicted 
from the gravity side at the leading order. 



Summary. — We have presented the first Monte 
Carlo calculations of the Wilson loop in a supersymmet- 
ric gauge theory at strong couphng. Up to subleading 
terms anticipated from the analysis on the gravity side, 
our results are in precise agreement with the prediction 
from the dual supergravity. This is a new and highly non- 
trivial evidence for the string/gauge duality. It would be 
nice to obtain the subleading terms explicitly from the 
gravity side, which will provide a nontrivial check of the 
duality including a' corrections. It is also interesting to 
extend this work to A/" = 4 SYM on R x S^, which is 



possible by using the equivalence 21| in the planar limit 
between the SYM and a mass-deformed MQM around a 
multi-fuzzy-sphere background. The equivalence is con- 
firmed by explicit calculations at weak coupling [23 |. 

The fact that we were able to see the Schwarzschild ra- 
dius of the dual black hole geometry by simulating large- 
N matrices gives us strong support and a firm ground for 
using matrix model simulations to study quantum grav- 



ity [23|. Note that the gauge theory description is valid 



also at small A and small N, where the dual supergravity 
description is no longer valid. Of particular interest is to 
study the parameter region corresponding to M-theory. 
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